Quasiparticle mass enhancement close to the quantum critical point in BaFe2(Asi_ x P x )2 
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We report a combined study of the specific heat and de Haas-van Alphen effect in the iron-pnictide super- 
conductor BaFe2(Asi_ a; P :E )2. Our data when combined with results for the magnetic penetration depth give 
compelling evidence for the existence of a quantum critical point (QCP) close to x — 0.30 which affects the 
majority of the Fermi surface by enhancing the quasiparticle mass. The results show that the sharp peak in the 
inverse superfluid density seen in this system results from a strong increase in the quasiparticle mass at the QCP. 
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Recently there has been much interest in the idea that super- 
conductivity in the iron-pnictides is boosted by the presence 
of a quantum critical point (QCP) located at the zero temper- 
ature end point of antiferromagnetic phase transition yj. It 
has been suggested that the presence of quantum critical fluc- 
tuations close to this point increase the superconducting tran- 
sition temperature by increasing the strength of the pairing 
interactions and / or by increasing the energy of the normal 
state relative to the superconducting state. 

In most iron-pnictides experimental evidence for quantum 
critical fluctuations has been difficult to find as its effects 
are easily masked by inhomogeneity and impurity scatter- 
ing. BaFe2(Asi_ :c P 2 ;)2 provides a particularly clean system 
in which to study these effects as the substitution of As by 
the isovalent ion P suppresses antiferromagnetism and induces 
superconductivity [2] without changing the electron/hole bal- 
ance and without inducing appreciable scattering [3,4]. These 
are likely the main reasons why the effect of the QCP on 
the physical properties are so clearly visible experimentally 
in BaFe2(Asi_ K P 2: )2 but not in its charge doped counterparts. 
Experimental evidence for quantum critical fluctuations in this 
system near the critical point x c ~ 0.3 include a funnel of 
T-linear behavior in the resistivity centered on x c |5[] and an 
increase in effective mass on one of the electron Fermi sur- 
faces, as measured by the de Haas-van Alphen (dHvA) effect 
IH, as x approaches x c . Nuclear magnetic resonance (NMR) 
experiments show that the magnetic ordering temperature ap- 
proaches zero at x c tffj. Most recently, measurements of the 
magnetic penetration depth A showed a large peak in A at x c 
indicating that the superfluid density is minimal at this critical 
x value and therefore that the QCP survives under the super- 
conducting dome |7|] . 

Although the above experiments have clearly shown that 
quantum critical fluctuations are present in this system close 
to x c , so far it is less clear if the quasiparticle mass renor- 
malization is the main factor which is causing the observed 
reduction in the superfluid density at x c fa] or whether or 
not it is significant on all Fermi surface sheets. Here we re- 
port a combined study of the specific heat and dHvA effect 
of BaFe2(Asi_ x P x )2 which indicates a strong increase in the 
effective mass in good quantitative agreement with the super- 
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FIG. 1: (color online) The inset shows the total measured specific 
heat for a sample of BaFe2(Asi-a;Pa;)2 with x — 0.30, the solid 
line is the fitted normal state background C„. The main part of the 
figure shows the specific heat with the normal state C n subtracted for 
different values of x and the solid lines show the construction used 
to determine the jump height AC and T c . 



fluid density measurements as suggested by recent theoretical 
work ]8[] . 

Single crystals of BaFe2(Asi_ a; P :E )2 were grown by a self- 
flux technique as described previously A particular re- 
quirement of both our specific heat and dHvA studies was to 
have very homogeneous samples. To this end we mostly fo- 
cused on using small crystals with masses 4-70 fig for specific 
heat and 0.1-4/ig for dHvA. To measure the specific heat of 
such small samples we used a custom build microcalorime- 
ter, where the thermometry elements and heater are deposited 
as thin films on a SiN membrane. Calibration of this device 
was checked by measuring high purity samples of Ag. We 
also used a more conventional long relaxation calorimeter |Ql 
for some larger samples (of order 300 /ig). In total around 30 
samples were measured using the microcalorimeter and only 
those with the very sharpest superconducting transitions were 
retained for further analysis. For these samples, we made 
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FIG. 2: (color online) The size of the jump in the specific heat 
AC/T C as a function of x. The line shows a fit to the logarithmic 
behavior expected near a QCP. The inset shows AC vs T c , with the 
solid line showing a fit to AC oc T™. The x dependence of T c is 
shown in the uppermost panel for the same samples and the dashed 
line is a guide to the eye. 



careful x-ray measurements to determine accurately the lat- 
tice constants and then used Vegard's law in conjunction with 
the lattice parameters for BaFe 2 As 2 and BaFe 2 P2 to deter- 
mine the value of x. dHvA measurements were performed 
by the torque method using micro piezoresistive cantilevers 
in both pulsed field (up to 60 T in Toulouse) and static field 
(up to 45 T in Tallahassee). Density functional theory (DFT) 
band-structure calculations were performed using Wien2K as 
described previously [10]. 

In Fig. Q] we show the jump in the specific heat C at T c for 
samples with x close to x c . Here the anomaly at the transition 
AC has been isolated from the phonon dominated background 
by subtracting a second order polynomial fitted above T c and 
extrapolated to lower temperature. Although there would be 
some uncertainty in using this procedure over an extended 
temperature range, the lack of appreciable thermal supercon- 
ducting fluctuations, as evidenced by the mean-field-like form 
of the anomaly, means that there is very little uncertainty in 
the size of AC. 

It is evident from the data in Fig. [T] that the size of the 
anomaly AC/T C depends very strongly on x and T c . This 
is shown in more detail in Fig. [2] As the strong increase in 



AC with x is accompanied by a relatively small increase in T c 
close to the critical point it seems highly unlikely that the cou- 
pling strength is the main factor. Indeed the normalized slope 
of C(T) - C n (T) just below T c , (T c /AC)d(C - C n )/dT, 
which depends strongly on the coupling streng th 111 ill , is al- 
most independent of x for 0.3 < x < 0.5 11211 . Instead 
we suggest that this increase in AC reflects the increase in 
the normal state density of states and hence the Sommerfeld 
coefficient 7. We show below that this interpretation gives 
remarkably good agreement with measurements of the quasi- 
particle mass from other techniques. We find that our data, for 
0.5 > x > 0.3, are well described by the logarithmic critical 
behavior, AC/T C = cq + c\ ln(x — x c ), expected close to a 

QCP d. 

A strong increase in AC with T c has been observed previ- 
ously in many different iron based superconductors including 
BaFe 2 (Asi_ x P a )2 HQ. Bud'ko, Ni and Canfield (BNC) 
II 1 3fl found that the data for materials with a wide range of T c 
could be described by the scaling law AC oc T c 3 . This has 
been interpreted as either originating from quantum critically 
1 1511 or from strong impurity pair breaking! 16] . In our sam- 
ples we find significant departures from this scaling and close 
to the critical doping a much stronger dependence on T c is ob- 
served; a fit to the power law AC oc T™ gives n = 6.5 ± 0.7 
for 30 K> T c >23 K (inset Fig.0. 

Although our absolute values of C (including phonons) are 
similar to those reported in Ref. II 1 411 . for samples with the 
highest T c we find values of AC which are up to a factor of 
two larger. The values of AC are very similar for the low 



3 



\ 2 - 

1 r- 



1 ( — I 1 1 1 - 

1.0-0.91 ln(x-0.30) 




o 




0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



FIG. 4: (color online) Dependence of the measured dHvA mass of 
the p orbits on x, the solid symbols are the new data from the present 
study and the open symbols are from our previous studies JH IToh . 
The solid line is a fit to the logarithmic behavior indicated. 



T c samples. This difference likely arises from the higher ho- 
mogeneity of the present samples. Indeed, we observed a 
similar reduction in AC for samples with broader transitions. 
The sharply peaked behavior of AC and strong asymmetry of 
T c (x) close to x c (Fig.O naturally leads to a strong suppres- 
sion of AC as the distribution of x within a sample becomes 
broader. 

Specific heat measures the total density of the states (DOS) 
and so contains contributions from all the Fermi surface 
sheets. According to density functional theory calculations 
ifioll for the end member BaFe2P2 the contributions from each 
sheet to the total DOS is 23% (band 1 hole), 37% (band 2 
hole), 20% (band 3 electron), 20% (band 4 electron). So the 
contribution from the holes is a little higher than from the 
electrons but they are roughly equal. dHvA effect measure- 
ments have the potential to resolve quasiparticle masses from 
individual orbits on each Fermi surface sheet. For BaFe2P2 
masses of almost all the observable orbits were reported by 
Arnold et al. ifioll . All orbits had relatively uniform mass en- 
hancements m* /mi, ranging from 1.6 to 1.9. For the As sub- 
stituted samples, the signals from the hole sheets are rapidly 
attenuated as x is decreased, but the signals from the inner (a) 
and outer (J3) electron sheets less so 13J. The j3 orbits are the 
most prominent and can be tracked to the lowest value of x 
and hence highest T c . In Fig. [3] we show how the amplitude 
of the dHvA signal from the /3 orbit reduces with temperature 
for several different values of x, along with fits to the standard 
Lifshitz-Kosevich (LK) formula 11711 which is used to deter- 
mine the orbitally averaged effective mass rri%. 

In Fig. |4] we show the effective mass of these orbits over 



FIG. 5: (color online) Mass enhancements as derived from 
specific heat, dHvA and microwave magnetic penetration depth 
measurements 1 7]. The dashed line is a guide to the eye. 



a wide range of doping, up to x = 0.38, T c — 28 K which is 
close to the maximum T c . For the lower values of x we are not 
able to resolve the maximal and minimal /3 orbits separately 
because of the restricted range of inverse field over which the 
signal is observable, so this mass represents an average of the 
two. The small size of these samples makes the accurate de- 
termination of the lattice constants very difficult. It was shown 
previously [3] that the frequency of both the a and (3 electron 
orbits decreases linearly with x, so we used this fact to deter- 
mine x. Because dHvA signals would be severely attenuated 
in any areas of the sample with inhomogeneous As/P content, 
this is a very accurate way to determine the exact x value for 
that part of the sample from which the dHvA signals originate. 
The solid line in Fig.|4]is a fit to the logarithmic form used for 
the specific heat. The point at x = 1 was not included in the 
fit as it is too far from the critical point. 

We now make a quantitative comparison between the dif- 
ferent measurements of effective mass. For the specific heat 
data, we first calculate 7 from the measured values of AC/T C , 
using the relation a — AC/jT c . Although a does vary with 
the strength of the electron-boson coupling and the anisotropy 
of the superconducting gap, it only increases by a factor of 
~2 even for extremely strongly coupled superconductors like 
Pb till] . Gap anisotropy decreases a, for example for d-wave 
a is 0.66 of the isotropic s-wave value. As both these effects 
are relatively small compared to the changes in AC(x) and 
will cancel each other out to some extent in BaFe2(Asi_ x P ;i; )2 
which has gap nodes [ 18], we make the approximation that a 
takes the weak coupling s-wave value of a = 1.43 for all x. 
Then we calculate an average effective mass enhancement by 
taking the ratio of 7 to the band-structure value calculated by 
DFT for BaFe 2 P 2 , j b = 6.94 mJ mol^K" 2 . 
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For the penetration depth data we make a similar compar- 
ison to the DFT calculations. Using the experimental values 
of Ao at low temperature measured by the microwave cavity 
method [7], we estimate the renormalization from m* x /mi, — 
AXl/Xl. Here A = n dHvA /n b = 0.358 + 0A4x is the mea- 
sured shrinkage of the Fermi surface volume with x 13] di- 
vided by the DFT volume for x = 1 and Af, = 660 A is 
the DFT(a; = 1) value of A. Although A is estimated from 
the electron sheet (a and f3) dHvA data only, the hole sheets 
must shrink similarly because the material remains compen- 
sated for all x. 

Finally, for the dHvA data we simply plot the measured 
values of m*o/mb for the j3 orbits. We use the (x = 1) DFT 
TOf, and jb values for all x for simplicity as these do not change 
appreciably with band energy shift corresponding to our range 
of x. 

In Fig. [5] we compare these estimations of the different av- 
erages of the effective mass. The data are remarkably con- 
sistent within experimental error over most of the phase di- 
agram. AC will be reduced at the edge of superconducting 
dome because of pair-breaking and sample inhomogeneity so 
we would not expect the samples with x = 0.2 and x = 0.65 
to provide a good estimate of m*. In principle, the mass en- 
hancements could vary significantly between different Fermi 
surface sheets. In the quasi-classical theory the superfluid 
density A~ 2 oc J v%.\v\ dS, and so the light electrons where 
the Fermi velocity v is high make the largest contribution. 
On the other hand, for the specific heat 7 oc J \v\~ 1 dS and 
so the heavy electrons contribute most. Hence the quantita- 
tive agreement between the trends for the mass enhancements 
indicated by A 2 , AC and the dHvA results for the electron 
/3 orbits tend to suggest that the mass enhancement is quite 
uniform over the Fermi surface. This does not exclude the 
presence of 'hot spots' or localized regions with much higher 
enhancement, which might result from regions that are par- 
ticularly well connected by strong spin-fluctuation modes, but 
does suggest rather similar average enhancement for the elec- 
tron and hole Fermi surfaces. This would be expected if the 
enhancement results from spin-fluctuated mediated scattering 
from the electron to hole surfaces with similar orbital charac- 
ter. These different experiments probe the quasiparticle mass 
in quite different regions of the (T, B) phase diagram. Pen- 
etration depth measures the mass at T = 0, B = 0, dHvA 
at T = but high B, and specific heat AC at B = and 
T = T c . The numerical consistency between AC/T C and A 2 
shown in Fig. therefore also suggests that the finite T of or- 
der T c does not significantly dampen the mass enhancements. 
The large B field used for the dHvA measurements also does 
not seem to have a strong effect except perhaps in the highest 
T c sample measurable by dHvA (x — 0.38). 

In summary, combined data for specific heat, de Haas-van 
Alphen effect and magnetic penetration depth [7] shows com- 
pelling evidence for the existence of a quantum critical point 
close to x = 0.30 in the BaFe2(Asi_ 2; P I )2 system, and that 
this effects the majority of the Fermi surface by enhancing 
the quasiparticle mass. The results show that the sharp peak 



in the inverse superfluid density seen in this system results 
from a strong enhancement of the quasiparticle mass at the 
QCP. The enhanced quasiparticle mass implies that the Fermi 
energy is suppressed, which is usually less advantageous for 
high T c . The fact that the highest T c is nevertheless attained 
right at x c = 0.3 with the most enhanced mass strongly sug- 
gests that the quantum critical fluctuations help to enhance 
superconductivity in this system. 
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